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I analyze the IR fixed point structure of Landau gauge QCD. Precisely the fixed point with a 
strong kinematic singularity of the quark-gluon vertex that proved crucial for the recently proposed 
confinement mechanism in the quenched approximation is absent in dynamical QCD. Therefore, 
the IR singularities do not induce asymptotic quark confinement but the long-range interaction is 
screened by unquenching loops at scales of the order of the quark mass. This provides the prerequisite 
for a microscopic description of deconfinement and string breaking. The fixed points determine the 
qualitative form of the heavy quark potential and may be relevant for hot and dense matter. 



Quark confinement is a manifestly non-perturbative 
phenomenon which is inherently scale dependent and fea- 
tures several crucial aspects that a possible mechanism 
has to explain. In order to give such an explicit dynam- 
ical confinement mechanism in terms of the underlying 
colored degrees of freedom this entails to actually specify 
them and correspondingly to fix a gauge. In the absence 
of dynamical quarks the most direct signature for confine- 
ment is an area law of large Wilson loops [I] that can be 
directly related to a linear rising potential between static 
color sources. In momentum space Green's functions this 
long-range physics is encoded in the infrared (IR) regime. 
The IR divergences of Green's functions are expected to 
offer an explanation for confinement since the early work 
of Weinberg [2]. In Landau gauge Mandelstam (3] pro- 
vided such a mechanism in a simplified approximation to 
the gluon Dyson-Schwinger equation (DSE) leading to a 
strongly divergent gluon propagator that scales ~ l/p 4 
in the IR limit. Yet, it has been shown recently that 
an IR enhanced gluon propagator is incompatible with 
the DSEs for the gluonic vertices |4J. Instead there is 
a consistent IR scaling fixed point of Yang-Mills theory 
(YMT) where the ghost sector of the theory dominates 
in the IR limit whereas the gluon dynamics is suppressed 
[5j |6]. This scaling solution provides a mechanism for 
gluon confinement within the scenarios of Kugo-Ojima 
[7] and Gribov-Zwanziger [8j[9]. ^ extends to arbitrary 
Green's functions [10J and kinematics |4J, represents the 
unique IR scaling fixed point [H HI] an d is also obtained 
with other functional methods However, it posed a 
puzzle how suppressed gluons can provide the strong in- 
teraction between color sources seen in lattice QCD [I]. 
Recently it was shown in the quenched approximation 
that this fixed point of Landau gauge YMT indeed in- 
duces a strong IR quark interaction that leads to a linear, 
confining potential between quarks [12]. This confine- 
ment mechanism relies on a self-consistently enhanced 
kinematic divergence of the quark-gluon vertex in the 
soft-gluon limit that the external gluon momentum van- 
ishes. The question is then how unquenching effects can 
circumvent this static confinement and thereby account 
for string breaking and hadronization seen in collider ex- 
periments. Such a mechanism must be related to scales 



of the order of the masses of created hadrons. In this 
work I outline a comprehensive analysis of the IR fixed 
point structure of Landau gauge QCD in the unquenched 
case. It confirms the screening of gluonic interactions via 
quark loops in the on-shell region that could "break the 
gluonic string" and drive the hadronization process. 
The non-perturbative dynamics of QCD is described by 
the system of DSEs which is generated algorithmically 
[13]. The DSEs form an infinitely coupled tower of equa- 
tions that in general requires approximations. Yet, it can 
be shown [14] that the equations for all n-point functions 
with n > 4 are linear in the sense that the correspond- 
ing n-point function appears at most once in each loop 
graph. As far as the IR singularities are concerned there 
is thereby no self-consistent enhancement mechanism in 
these equations. They also include no tree level term and 
are hence effectively determined by lower order Green's 
functions. Therefore, any IR enhancement has to be gen- 
erated by the primitively divergent Green's functions and 
the truncation can be restricted to the latter. Moreover, 
it has been found previously |4, 10J that all contributions 
including a 4-gluon vertex are strongly IR suppressed 
so that this vertex can be neglected as well. However, 
the DSE system features a peculiarity compared to other 
functional approaches like the functional renormalization 
group (FRG) [llj or n-particle irreducible (nPI) actions 
[T5] in that it involves a bare vertex in every DSE so that 
IR strength can be missing in the lowest order diagrams 
when they involve enhanced vertices. Since ghosts and 
quarks do not couple at the tree level this requires to 
include the ghost-quark vertex to consider all IR leading 
terms. The considered truncation is then given by the 
DSEs in the gauge and matter sector in figs. [I] and [2] 
For the analysis of the IR regime a flexible power count- 
ing scheme [4j [6j [10] has been developed that yields the 
IR fixed point structure without an actual solution of 
the non-linear integral equations. Here I only sketch the 
main steps of the analysis which is similar to the one pre- 
viously performed in the gauge sector |4J and extends the 
one in the quenched approximation [12]. The details will 
be presented in a forthcoming article [14]. The infrared 
analysis relies on the observation that YMT is scale in- 
variant at the classical level. Quantum fluctuations in- 
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Figure 1: The gauge part of the considered coupled DSE sys- 
tem. Ellipses denote graphs involving neglected 4- and 5-point 
functions and closed quark loops are absent in quenched QCD. 




Figure 2: The matter part of the considered DSE system. 

duce a dynamical scale Aqcd and in QCD the matter 
sector introduces in addition explicit mass scales, but it 
is experimentally established that strong interaction has 
no scales far below the MeV scale. Correspondingly, it is 
a reasonable assumption that in the IR regime p<C Aqcd 
the scale dependence of all Green's functions is given by a 
power law scaling in terms of the the dimensionless ratio 
P=p/Aqcd, parametrized for the propagators by 
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where S g i, S g h and S q denote the anomalous IR exponents, 
distinguished by lower indices. The vertices depend on 
several external momenta and can feature divergences in 
various kinematic configurations that are described by 
different IR exponents @J[T2], distinguished by upper in- 
dices. In the uniform limit that all external momenta 
vanish these exponents are S gg for the ghost-gluon ver- 
tex, (Jg f° r the 3-gluon vertex, 5^ g for the quark-gluon 
vertex and 5™ h for the quark-ghost vertex. In addition 
there are exponents dqgh wnen two external 

momenta vanish, 5$, 5* 5$* 5% Q1 5 9 Q \, S q Qah when a 



single external momentum vanishes and finally 8 s qgh when 
only the difference of the two external quark respectively 
ghost legs is soft. Where there are different tensor struc- 
tures these are the exponents of the IR leading part(s). 
The loop integrals in the DSEs are dominated by the sin- 
gularities of the integrand and receive contributions from 
loop momenta in the vicinity of all external momentum 
and mass scales [4], i.e. in QCD also by scales of the order 
of the constituent quark masses. In the IR limit there is 



a clear scale separation and one can divide the contribu- 
tions from soft and hard loop momenta. By determining 
the IR scaling of these distinct contributions from the 
different loop graphs via a power counting analysis the 
DSE system translates to a coupled set of algebraic equa- 
tions for the anomalous exponents. In each DSE the IR 
leading term on the right hand side, corresponding to the 
smallest IR exponent, determines the IR behavior of the 
Green's function. When hard momenta are present there 
are various possible cancellations in the IR analysis that 
affect the canonical scaling of Green's functions. These 
arise from gluon transversality, the symmetry of Green's 
functions [4] and from the Dirac structure [12J. There is 
an additional suppression of the contribution from hard 
loop momenta k where the quark propagator can be ex- 
panded in the soft external momentum p 
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which shows that if there is no constituent quark mass 
m q = M(0) = or another hard scale in the integral the 
leading term is scale-free and is removed in the renor- 
malization procedure so that the propagator is actually 
suppressed in p. The gauge propagators feature an analo- 
gous suppression by p 2 [4]. In the power counting analysis 
this is considered by symbols that take values H q = \ and 
figi = 1 when the propagators are massless, i.e. 5 q = — - 
and S g i =0, and which vanish otherwise. Analogous sym- 
bols mP q and jj/q are introduced at the tree level. 
The power counting analysis for the DSE system figs. 
[I] and [2] results in an extensive system of coupled al- 
gebraic equations for the above 18 IR exponents [14]. 
Fortunately, there are important constraints on the IR 
exponents, given in [4], and analogous ones in the quark 
sector that simplify the system considerably. Previously 
such constraints had been obtained from the assumption 
that a skeleton expansion of the vertices should not ex- 
plicitly diverge. Yet, the same constraints are obtained 
independently from the complementary set of FRG equa- 
tions without additional assumptions [TTJ ITS). Since the 
equations for the quark-ghost vertex are linear they can 
be solved by successive replacements of the possible so- 
lutions and the connection between the exponents in dif- 
ferent kinematic limits provides further complementary 
constraints [14]. Applying the constraints yields that the 
fermionic 3-point vertices are not divergent when only a 
fermion momentum vanishes 5^ g = 5% g = and leaves the 
simplified system for the IR exponents of the propagators 
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and of the remaining primitively divergent vertices 
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where the contributions from different graphs in figs, [I] 
and [2] are separated by semicolons and contributions from 
different regions of the same loop integral by commas. It 
is instructive to recall the solution of the quenched sys- 
tem [12] first, where the terms in brackets arising from 
closed quark loops are absent. There the quark-gluon ver- 
tex features different solutions with a trivial respectively 
a strongly enhanced IR scaling. In the DSEs the latter 
arises from the IR region of the ghost loop diagram in fig. 
[2] corresponding to the first non-trivial term in eqs. (|8| 
and ([9]). Inserting the dressed quark-ghost vertex, results 
in a non-Abelian graph with a ghost-triangle, represent- 
ing the IR leading contribution |4J to the 3-gluon vertex 
in eq. ([6|. Since such a fully dressed non-Abelian graph 
arises directly both in the FRG [11] and nPI actions [15], 
the dynamics seems to be represented more efficiently in 
these approaches. More importantly, for the strong soft- 
gluon singularity the unquenching contributions in the 
last term in brackets in eqs. (|3| and (|6| are incompatible 
with the gauge sector, e.g. 1 — S g i ^ 5|^. 
The starting point for the full solution of the unquenched 
system is eq. ([I]) for the ghost propagator. As it stands 
the bare term is leading and in addition to the trivial, 
perturbative solution that is also realized in the ultravi- 
olet (UV) regime one obtains a decoupling solution with 
a bare ghost and an IR finite gluon propagator [TTJ [18]. 
Due to the suppressed gluonic interaction there are no 
IR enhanced Green's functions. The decoupling solution 
seems to be found in current lattice studies [19]. Yet 
the study of the IR regime of a gauge theory on the lat- 
tice is a notoriously hard problem and there are several 
issues like the present implementation of the gauge con- 
dition and the related contamination with Gribov copies 
that demand caution [18J. In particular, since there is 
another solution that explains important qualitative as- 
pects of strong interaction physics better: For an IR 
enhanced propagator the renormalization procedure can 
also be performed in a way that the tree level term in 




Figure 3: left: The DSE for the 4-quark vertex and IR lead- 
ing contributions to the arising 4-quark-gluon vertex, right: 
Heavy quark potential with areas of application of the pertur- 
bative (V) UV, static (<S) and massive (M) IR fixed points. 



the ghost equation is cancelled identically. This yields 
the scaling solution where S g i = —25 g h = 2n |5j, depend- 
ing on a parameter n > since the gluon equation (|3| 
is then identically fulfilled by the ghost-loop exponent 
[4]. The correspondingly subleading contribution from 
unquenching loops with hard momenta in eq. ^ yields 
the important constraint 5* > 1— 2k— fi q . The scaling so- 
lution depends on the three qualitatively different cases 
of an infinite, a finite or a vanishing constituent quark 
mass m eq. The first case is similar to the quenched 
limit, but here quarks are entirely static and merely act 
as sources. Using the available constraints the solution of 
the residual system is possible without further assump- 
tions and the resulting IR fixed points of QCD are given 
in table [T| Strikingly for all QCD fixed points the gauge 
sector is not altered by the quark dynamics compared to 
the YMT result [1[10]. The main result of this article is, 
however, that due to the above constraint the solution 
with a strong soft-gluon singularity of the quark-gluon 
vertex is ruled out in dynamical QCD and only the so- 
lution with a trivial soft-scaling remains. This happens 
since hard virtual quark loops with momenta around the 
quark mass screen the long-range gluonic interaction. 

Now let us discuss the impact of the IR divergences on 
quark confinement which should in the case of mesons 
be described by the 4-quark vertex DSE given in the first 
line of fig. [3J The quarks in hadrons have momenta of the 
order of hundreds of MeV and correspondingly the rele- 
vant kinematics for confinement is given by hard external 
quark momenta. Though the exchanged gluon momen- 
tum becomes soft when they are far spatially separated. 
Using the IR exponents in the static limit given in ta- 
ble [TJ the first diagrams in fig. [3] do not induce strong 
long-range interactions since IR strength is missing due 
to the bare quark-gluon vertex [12]. However, this is not 
the case for the last graph with the higher order 4-quark- 
gluon vertex. Its DSE in the second line of fig. [3] includes 
graphs where all four quark-gluon vertices attached to the 
exchanged soft gluons are dressed. The kinematic regime 
where the small loop in the first line is dominated by mo- 
menta of the order of the quark mass and the one in the 
second line is soft yields a far stronger divergence. In the 
heavy mass limit the hard loop shrinks to a point and the 
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Table I: The anomalous power law exponents of the leading Green's functions for the different fixed points of QCD distinguished 
by the constituent quark mass m q . These power laws are only valid up to logarithmic corrections. The full vertices include 
also the canonical scaling dimension —1 for the gauge propagators and - for the gauge 3-point vertices in the uniform limit. 
The value of k is fixed by an explicit IR solution and the best known value is ft ~ 0.59 [6J. There are additional fixed points in 
the chiral case where the quark propagator is IR divergent [TJ, yet physically they are excluded by the current quark masses. 



leading corrections reduce to graphs similar to the first 
two in the 4-quark DSE but with all four quark- gluon 
vertices dressed - again in complete agreement with the 
FRG [II] or nPI methods [15]. With the strong soft- 
gluon singularity of the dressed quark-gluon vertices this 
yields an IR interaction that scales ~l/p 4 corresponding 
to a linear, confining potential in coordinate space [12J. 
In contrast, in the unquenched case the quark-gluon ver- 
tex is not enhanced corresponding to a 4-quark interac- 
tion ~p 8K which yields a strongly suppressed long-range 
behavior and cannot be interpreted by a non-relativistic 
potential. Instead, the long-range interaction, found in 
the static respectively quenched case, is screened by vir- 
tual quark loops at scales of the order of the quark mass 
and dynamical QCD is indeed not asymptotically con- 
fining. Yet, in the heavy mass limit m q ^> Aqcd and 
for momenta p<^m q where the quenched limit is a valid 
approximation one can expect a linear potential over dis- 
tances for which V ~ ax <C 2m q . In the short distance 
regime p~1/x^>Aqcd the perturbative UV fixed point 
is realized which yields an interaction ~l/p 2 correspond- 
ing to a Coulomb potential, whereas the potential picture 
breaks down for distances x < l/m q . Correspondingly 
the qualitative form of the heavy quark potential is en- 
tirely determined by the fixed points of QCD as shown 
schematically in fig. [3j The arising screening scale is pre- 
cisely of the order of the masses of corresponding mesonic 
states and thereby it is not hard to imagine that "the 
string breaks" and real particles are produced when the 
system is sufficiently excited. This is described by higher 
Green's functions, like the 6-quark vertex, and it remains 
to be shown that only color singlets are produced. How- 
ever, due to the intermediate rise in fig. [3] it is clear that 
produced quarks are confined into bound states when 
their excitation energy is small enough. The presence of 
light quarks surely goes beyond the present IR analysis 
and requires computations at finite hadronic scales. 
Finally the above fixed points could be relevant for hot 
and dense matter. In the limit of asymptotic temperature 
T and/or chemical potential \i the quark masses become 
irrelevant in comparison and the chiral cases in table [T| 
which do not induce quark confinement, may be relevant 



for some scale region. Yet, in addition to the perturba- 
tive fixed point, describing complete deconfinement, the 
chiral scaling fixed points feature a confined gauge sector. 
The stronger one exhibits a Coulombic quark interaction 
as in the perturbative case and may be relevant for the 
strongly coupled plasma observed in heavy ion collisions. 
In contrast, the weaker one corresponds to a quark sector 
that completely decouples from a strongly coupled gauge 
sector in the IR limit. Interestingly, such a behavior has 
been independently found in an analysis of the low energy 
dynamics of ungapped modes at high \i |2Q| . 
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